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Abstract 

The van Vleck determinant is an ubiquitous object, arising in many 
physically interesting situations such as: (1) WKB approximations to 
quantum time evolution operators and Green functions. (2) Adiabatic ap- 
proximations to heat kernels. (3) One loop approximations to functional 
integrals. (4) The theory of caustics in geometrical optics and ultrason- 
ics. (5) The focussing and defocussing of geodesic flows in Riemannian 
manifolds. 

While all of these topics are interrelated, the present paper is par- 
ticularly concerned with the last case and presents extensive theoretical 
developments that aid in the computation of the van Vleck determinant 
associated with geodesic flows in Lorentzian spacetimes. A fortiori these 
developments have important implications for the entire array of topics 
indicated. 
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1 INTRODUCTION 



The van Vleck determinant is a truly ubiquitous object in mathematical physics. 
Introduced by van Vleck in 1928 0, it was first utilized in elucidating the 
nature of the classical limit of quantum mechanics via WKB techniques. Further 
developments were due to Morette (2| . A nice discussion of this original line of 
development can be found in Pauli's lecture notes ||. 

In a different vein, de Witt j|, ||, |(| [t]] and others have developed extensive 
and powerful adiabatic expansion methods for approximating quantum Green 
functions in Lorentzian and Ricmannian manifolds. Textbook discussions may 
be found in Birrell and Davies || and Fulling || . 

By extension, these adiabatic techniques lead to powerful point -splitting 
techniques. These techniques are useful, for example, for estimating the vac- 
uum expectation value of the renormalized stress-energy tensor of a quantum 
field propagating on a Lorentzian spacetime. As will soon be apparent, in this 
particular context the van Vleck determinant is essentially a measure of the 
tidal focussing and/or defocussing of geodesic flows in spacetime. 

This paper will develop several formal techniques for getting a handle on 
the van Vleck determinant. The evolution of the van Vleck determinant as one 
moves along a geodesic will be explicitly evaluated in terms of the expansion 
of the geodesic flow. The relationship with the (spacelike/lightlike/timelike) 
versions of the Raychaudhuri equation is thus made manifest. Constraints, 
such as the Weak Energy Condition (WEC), are then used to develop powerful 
inequalities limiting the behaviour of the van Vleck determinant. The evolution 
of the van Vleck determinant is then reformulated in terms of tidal effects — the 
evolution of a set of Jacobi fields under the influence of the full Riemann tensor. 
A formal solution to the resulting integral equation is presented. Using these 
formal techniques, a weak field approximation is developed. This weak field 
approximation will be seen to depend only on the Ricci tensor of spacetime. 
In a similar vein, a short distance approximation is developed. Finally, an 
asymptotic estimate is provided of the behaviour of the determinant as one 
moves out to infinity. 

The physical problem that originally stimulated my interest in these inves- 
tigations was Hawking's Chronology Protection Conjecture |h], pl| . That con- 
jecture led to several different computations of the vacuum expectation value of 
the renormalized stress-energy tensor for quantum fields propagating in space- 
times on the verge of violating chronology protection (lj, [l3| |l4|, [l5| . The van 
Vleck determinant is an overall prefactor occurring in all those computations. 

It is nevertheless true that much of the mathematical machinery developed 
in this paper has a considerably wider arena of applicability. Accordingly, some 
effort will be made to keep the discussion as general as reasonably possible. 

Notation: Adopt units where c = 1, but all other quantities retain their 
usual dimensionalities. In particular G = h/m 2 p = i^/fr. The metric signature 
is taken to be (— , +,•••, +). General conventions follow Misner, Thorne, and 
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Wheeler @. 



2 THE VAN VLECK DETERMINANT 

2.1 General definition 

Consider an arbitrary mechanical system, of n degrees of freedom, governed by 
a Lagrangian C(q,q). Solve the equations of motion to find the (possibly not 
unique) path 7 passing through the points (qi,U) and (qf,tf). Then calculate 
the action of that path 

S-y{qi,U;qf,t f ) = C(q,q)dt. (1) 
The van Vleck determinant is then defined as 

/Mft, W/ ) 5 (-1)" det { d2S ^XT tf) } (2) 

The van Vleck determinant, in this original incarnation, occurs as a prefactor 
in the WKB approximation to the quantum time evolution operator |Q 
In the Schrodinger picture 

< qf\exp{—iH(tf ~ U)/h}\qi >« 

{2mh)~ n l 2 \l A 7fe, Ul qf, tf) exp{+iS 7 ( qi , U1 qf, tf)/h}- (3) 

7 

In the particular case of a geodesic flow on a Lorentzian or Riemannian manifold 
the appropriate action to be inserted in this definition is the geodetic interval 
between a pair of points. 

2.2 Geodetic interval 

Consider a Lorentzian spacetime of total dimensionality (d + 1). That is: d 
space dimensions, 1 time dimension. With appropriate changes, these results 
apply also to a Riemannian manifold of total dimensionality (d+l). To go from 
the Lorentzian to the Riemannian case the only significant change is that one 
should forget all the subtleties associated with lightlike (null) geodesies. 
The geodetic interval may be defined by: 

a 1 (x,y) = ±^[s 1 {x 1 y)} 2 (4) 

Here we take the upper (+) sign if the geodesic 7 from the point x to the point 
y is spacelike. We take the lower (— ) sign if this geodesic 7 is timelike. In either 
case we define the geodesic distance s 1 (x,y) to be positive semi-definite. 
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Note that, provided the geodesic from x to y is not lightlike, 



V£cr 7 (x,y) = ±s J (x,y) V*s 7 (x,y), (5) 
= +s 7 (x, y) t^(x; 7; x <- y), (6) 

= ( 7 ) 

Here i M (x;7;x <— y) = ±g A " y VJ;s 7 (a;, y) denotes the unit tangent vector at the 
point ir pointing along the geodesic 7 away from the point y. When no confusion 
results we may abbreviate this by t^(x <— y) or even t^{x). 

If the geodesic from x to y is lightlike things are somewhat messier. One 
easily sees that for lightlike geodesies V^a 7 (x,y) is a null vector. To proceed 
further one must introduce a canonical observer, described by a unit timelike 
vector V* 1 at the point x. By parallel transporting this canonical observer 
along the geodesic one can set up a canonical frame that picks out a particular 
canonical affine parameter: 

V^o- 7 (a;, y) = +C 7 0, y) l?(x; 7; x <- y), (8) 
l ll (x;r,x^y)V' i = -l, (9) 
Q 1 {x,y) = -V^a 1 {x,y). (10) 

Note that this affine parameter £ can, crudely, be thought of as a distance 
along the null geodesic as measured by an observer with four-velocity V^. By 
combining the null vector I 11 with the timelike vector V 1 , one can construct a 
second canonical null vector: = 2V^ — I 11 . A. quick calculation shows 

V*C 7 0,y) = ~\ m^{x;r,x <- y), (11) 
m„(i) V = -1; (12) 
m^ar) = -2. (13) 

Unfortunately, while spacelike and timelike geodesies can be treated in a unified 
formalism as alluded to above, the subtleties involved with lightlike geodesies 
will require the presentation of several tedious variations on the general analysis. 

2.3 Specific definition: geodesic flow 

Consider geodesic flow in a Lorentzian spacetime. In the present context the 
(scalarized) van Vleck determinant is defined by 

V9Wg(y) 

This definition has a nice interpretation in terms of the Jacobian associated with 
the change of variables from (x,t) to (x,y). One may specify a geodesic either 
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by (a) specifying a single point on the geodesic, x, and the tangent vector, t, 
at that point; or by (b) specifying two separate points, (x,y), on the geodesic. 
The (scalarized) Jacobian associated with this change of variables is 



j {x = det{d(x,t)/d(x,y)} = dee\dt/dy} 
^9{x)g{y) \/g{x)g{y) 



Here the det' indicates the fact that we should ignore the known trivial zero 
eigenvalue(s) in determining this Jacobian. 

If the points x and y are spacelike or timelike separated then the one trivial 
zero arises from the fact that the tangent vector is normalized, t^t^ = ±1, so 
that t^(x)\7f,t^(x) = 0. To sec the connection with the van Vleck determinant, 
observe 

V*V> 7 (x,y) = VZ[s 1 (x,y)t ll (x;r,x*-y)], 

= 8^(x,y)Vlt li (x;r,x<-y)±t v (y)t li (x). (16) 

By adopting suitable coordinates at x, and independently at y, one can make the 
tangent vectors t M (x) and t M (y) both lie in the "t" direction (if they are timelikc) 
or both lie in the "z" direction (if they are spacelikc). Then by inspection 

A l{ x, y) = ±[-< h ix,v)]< ^Y;:Tr y)} 

V9{x)g{y) 

= ±l-s^x,y)] d J 7 (x,y). (17) 

If the points x and y are lightlike separated then the two trivial zeros arise 
from the fact that the tangent vector satisfies two constraints: l^l^ = 0, 
and Z*% = -1. Thus both l» {x)V v v t ^(x) = 0, and V^{x)W v v t^{x) = 0. To 
precisely determine the connection with the van Vleck determinant requires 
some subtlety For a point x close to, but not quite on, a null geodesic emanating 
from y one may usefully decompose the gradient of the geodetic interval as 

V x ^i(x,y) = Hj(x,y) lp(x;r,x <- y)+^(x,y) m M (a;;7;x <- y). (18) 

Here ( and £ are to be thought of as curvilinear null coordinates. They are the 
curved space generalizations of (t ± x)/2. One is now in a position to calculate 

V*V»o- 7 (x,y) = V v l/ [((x,y)l t ,(x^y)+((x,y)m f ,(x^y)}, 

= C(x, y) V^(a;; 7; x <- y) + £(x, y) V^m M (x; 7; x <- y) 

-^ m Ay) iA x ) - 2 l »(y) ( 19 ) 

Now go to the light cone, by setting £ = 0. By adopting suitable coordinates at 
x one can arrange: l^{x) = (1, 1, 0, 0) and m fl (x) = (1, —1, 0, 0). Independently, 



5 



one can arrange the same to be true at y. Finally, careful inspection of the 
above reveals 

V9{x)g{y) 

= (-l)%(x,y)]V-» J^x,y). (20) 

2.4 Elementary results 

In view of the equation 

V^cr 7 (a;,y) = s 7 (x,y) t^x-^^x <- y), (21) 
one has as an exact result (valid also on the light cone) 

a 7 (x,y) = -g"" V*a 7 (x,y) V> 7 (z,y). (22) 

Repeated differentiations and contractions result in j|, ||, || [?J 

V£ [A 7 (x, y) V*a 7 (x, y)] = (d+ l)A 7 (x, y). (23) 

Assuming that x and y are either timelike separated or spacelike separated, this 
may be rewritten as follows 

V£ [A 7 (z, y) s y (x, y) t^x)] = (d + 1)A 7 (.t, y). (24) 

Use of the Leibnitz rule leads to 

s 7 (x,y) t^(x) V£A 7 (x,y) 
±t ll (x) i?{x) A y (x,y) 
+s 7 (x,y) A 7 (x,y) V x %(:r) 

= (d+l)A 7 (z,y). (25) 

One notes that defines a normalised spacelike or timelike vector field. 

Thus 

t^x) V£ A 7 (x, y) = (^zj^y ~ (V/t,(x))j A 7 (x, y). (26) 

Now define A 7 (s) to be the van Vleck determinant calculated at a proper dis- 
tance s along the geodesic 7 in the direction away from y and towards x. Note 
the boundary condition that A 7 (s = 0) = 1; and define 

t»( x )Vlf = f s . (27) 
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One finally obtains a first order differential equation governing the evolution of 
the van Vleck determinant 



dA 7 (s) (d 
ds \ s 



A 7 ( s ). (28) 



Here is the expansion of the geodesic spray defined by the integral curves of 
t^{x) : that is, 9 = V^t^. Direct integration yields 



A 7 (a:,y) = s 1 (x, y) d exp ^~ 
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6ds (29) 



Here the integration is to be taken along the geodesic 7 from the point y to the 
point x. The interpretation of this result is straightforward: take a geodesic 
spray of trajectories emanating from the point y. If they were propagating in 
flat space, then after a proper distance s the (relative) transverse density of 
trajectories would have fallen to s~ d . Since they are not propagating in flat 
spacetime the actual (relative) transverse density of trajectories is given by the 
exponential of minus the integrated expansion. The van Vleck determinant is 
then the ratio between the actual density of trajectories and the anticipated flat 
space result. Note that by explicit calculation, the van Vleck determinant is 
completely symmetric in x and y. 

An analogous development holds if the points x and y are connected by a 
null geodesic. The details are, as one has by now grown to expect, somewhat 
tedious. One returns to equation d23| ) and uses equation (|l^) to derive 

V^A 1 ( I ,j / ){(/ (I (a;)+(m,(i)}]=(d+l)A 7 (a;, ! /). (30) 

Use of the Leibnitz rule, followed by the limit £ — * 0, now leads to 

CJ M (aO V£A 7 (x,2/) 
-l^x) m^ix) A y (x,y) 
+C A 7 (x,y) V x %{x) 

= (d + l)A 7 (x,y). (31) 

One recalls that, by definition, l^{x)m^{x) = —2. Adopting suitable definitions, 
analogous to the non-null case, one derives 



dA 7 (C) f[d-l] 



dC V c 

Direct integration yields 



A 7 (C). (32) 



A 7 (x, y) = C 7 (ar, y) [d ~ 1] exp (- j Od^j (33) 
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In particular, note that while the affine parameter £, and the expansion 9 = 
V^/^ both depend on the canonical observer V 11 , the overall combination is 
independent of this choice. The appearance of the exponent d — 1 for null 
geodesies, as contrasted to the exponent d for spacelike or timelike geodesies 
might at first be somewhat surprising. The ultimate reason for this is that in 
(d + 1) dimensional spacetime the set of null geodesies emanating from a point 
x sweeps out a d dimensional submanifold. As one moves away from x an affine 
distance ( into this submanifold the relative density of null geodesies falls as 
C 1 ■ In (3 + 1) dimensional spacetime, this is just the inverse square law for 
luminosity. 

Equation (^2|) above should be compared with equation (22) of Kim and 
Thorne The typographical error in that equation (the 3 should be a 2) 

fortunately does not propagate into the rest of their paper. The virtue of this 
tedious but elementary analysis is that one now has a unified framework appli- 
cable in all generality to cases of timelike, lightlike, or spacelike separated points 
x and y. 



2.5 Inequalities 

Now that one has these general formulae for the van Vleck determinant, powerful 
inequalities can be derived by applying the Raychaudhuri equation and imposing 
suitable convergence conditions. Restricting attention to vorticity free (spacelike 
or timelike) geodesic flows the Raychaudhuri equation particularizes to 

^ = -(i^ r) _ 2CT 2 * (34) 
ds d 

For a concrete textbook reference, this is a special case of equation (4.26) of 
Hawking and Ellis |l7j] , generalized to arbitrary dimensionality. While Hawking 
and Ellis are discussing timclike geodesies, the choice of notation in this paper 
guarantees that the discussion can be carried over to spacelike geodesies without 
alteration. The quantity a denotes the shear of the geodesic congruence, and 
a 1 is guaranteed to be positive semi-definite. 

Definition: A Lorentzian spacetime is said to satisfy the (timelike/null/spacelike) 
convergence condition if for all (timelike/null/spacelike) vectors t M : [R^t^t") > 
0. ■ 

Definition: A Riemannian manifold is said to possess semi-positive Ricci curva- 
ture if for all vectors t^: (R^W) > 0. ■ 

If the points x and y are (timelike/spacelike) separated, and the spacetime 
satisfies the (timelike/spacelike) convergence condition then one has the inequal- 
ity 

dd e 2 

Ts * ~~d' (35) 
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This inequality is immediately integrable 



6^)-W) + S d- (36) 

For the geodesic spray under consideration, one has 0(s) — > d/s as s — > 0. 
Therefore [1/0(0)] = 0. One has thus derived an inequality, valid for geodesic 
sprays, under the assumed convergence conditions 

6(a) < +~. (37) 

This implies that the van Vleck determinant is a monotone function of arc 
length, dA/ds > 0, and immediately integrates to an inequality on the deter- 
minant itself A 7 (x, y) > 1. 

For null geodesic flows the Raychaudhuri equation becomes 

| = - 2 " 2 -(^T)' < 38 » 

See equation (4.35) of Hawking and Ellis |l7j]. If the spacetime now satisfies the 
null convergence condition one can deduce the inequality 

— > . (39) 

d( ~ d-1 v ' 

Upon integration, and use of the relevant boundary condition, one derives in a 
straightforward manner the inequality 

0(0 < +^r^- (40) 

Which implies that the van Vleck determinant is a monotone function of the 
null affine parameter, g?A/g?£ > 0, and immediately integrates to an inequality 
on the determinant itself A 7 (x, y) > 1. One is now in a position to enunciate 
the following results. 

Theorem: In any Lorentzian spacetime, if the points x and y are (timelike/null/ 
spacelike) separated, and the spacetime satisfies the (timelike/null/spacelike) 
convergence condition, then the van Vleck determinant is a monotone function 
of affine parameter and is bounded from below: A 7 (x, y) > 1. ■ 
Theorem: In any Riemannian manifold of semi-positive Ricci curvature, the van 
Vleck determinant is a monotone function of affine parameter and is bounded 
from below: A 7 (x,y) > 1. ■ 

To see the physical import of the (timelike/null/spacelike) convergence con- 
ditions, consider a type I stress-energy tensor. The cosmological constant, if 
present, is taken to be subsumed into the definition of the stress-energy. Then 
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Tfj, v ~ diag[p;pi]^i/, where the pi are the <i principal pressures. From the Ein- 
stein equations, i? A ,„ — \Rg^v — 8ttG T^ v , one deduces R^ v ~ 47rG diag[p + 

E i ^;P + 2 K -Ej Pjl^- 
The timelike convergence condition implies: (i) p+J2jPj — 0; (2) Vi,p+Pi > 
0. This is equivalent to the Strong Energy Condition (SEC). 

The null convergence condition implies: Vi,p + Pi > 0. This is a somewhat 
weaker restriction than the Weak Energy Condition (WEC). (The WEC would 
require the additional constraint p > 0.) 

The spacelike convergence condition implies: (1) Mi, p + 2p t — J2j Pj > 0; (2) 
\/i,p + Pi > 0. This is a rather strong constraint that is not equivalent to any 
of the standard energy conditions. Particularize to (3+1) dimensions, then the 
constraints (1) can be made more explicit as: p+pi— pi— pz > 0; p+pi—pz—p\ > 
0; p + P3 — pi — P2 > 0. By adding the first two of these constraints one may 
deduce p — P3 > 0. Similarly Vi,p — pi > 0. Thus the spacelike convergence 
condition implies the Dominant Energy Condition (DEC), which in turn implies 
the Weak Energy Condition (WEC) . The implication does not hold in the reverse 
direction. In particular, for an electromagnetic field of zero Poynting flux T^ v ~ 
p ■ diag[+l; — 1; +1; +1]^. This stress-energy tensor satisfies the Dominant 
Energy Condition but not the spacelike convergence condition. A canonical 
example of a stress-energy tensor that does satisfy the spacelike convergence 
condition is a perfect fluid that satisfies the Dominant Energy Condition: p+p > 
0; p -p > 0. 

Note that the timelike and spacelike convergence conditions separately imply 
the null convergence condition. 

The physical reason behind the inequality satisfied by the van Vleck determi- 
nant is now manifest: "Ordinary" matter produces an attractive gravitational 
field. An attractive gravitational field focuses geodesies, so that they do not 
spread out as much as they would in flat space. Then the van Vleck determi- 
nant is bounded from below, and in fact continues to grow as one moves along 
any geodesic. 

One should note however, that several interesting classes of spacetimes vio- 
late one or more of these convergence conditions. This makes the inequalities 
derived here somewhat less useful than they might otherwise be, and suggests 
that it would be useful to explore the possibility of deriving more general inequal- 
ities based on weaker energy conditions. For example, spacetimes containing the 
interesting class of objects known as traversable wormholes violate all three of 
these convergence conditions |L8[ [l^]. Indeed in that particular class of space- 
times the required presence of "exotic" matter leads to van Vleck determinants 
that are arbitrarily close to zero jl2| . 

2.6 Reformulation: Tidal focussing 

Direct computations using the preceding results are unfortunately rather dif- 
ficult. To get a better handle on the van Vleck determinant we shall first 



10 



apparently make the problem more complicated by deriving a second order dif- 
ferential equation for the van Vleck determinant. By going to the second order 
formalism it is possible to relate the van Vleck determinant directly to the tidal 
forces induced by the full Riemann tensor — more specifically to the focussing 
and defocussing effects induced by tidal forces. 

The most direct route to this end is to pick up a standard reference such as 
Hawking and Ellis |l7). Equation (4.20) on page 83 implies that 

J 6ds = bxdet[Af i v (s)]. (41) 

Here A^ v {s) is the d x d matrix describing the evolution of the separation of 
infinitcsimally nearby geodesies. This formalism may also be reformulated in 
terms of the Jacobi fields associated with the geodesic 7. Note that the d x d 
matrix A fJ - l/ (s) may equivalently be thought of as a (d+ 1) x (d+ 1) matrix that 
is trivial in the extra entries: Au+i) — Ad® s. In terms of A 

A 7 {x, y) = s 7 (z, y) d det^ 1 ) = det(s 7 (a;, y)A~ 1 ). (42) 

Where the last determinant can be taken in the sense of either a d X d or a 
(d + 1) X (d + 1) matrix. 

An alternative, more explicit, but also more tedious route to establishing 
the preceding equation is to define the object 

A(s) ee Pexp{ /(V & t)ds}. (43) 

J "V 

Here the symbol P denotes the path ordering process. By definition of path 
ordering 

dA -> f 

— = (V <g> t) ■ Pexp{- / (V ® t)ds}. (44) 
ds J 7 

On the other hand, A(0) = I. Comparison with equation (4.10) of Hawking and 
Ellis now shows that A(s) = A(s) tl v <g> dx v . Finally note that 

det A = exp{ / tr(V ® t)ds} = exp{ / 6ds}. (45) 
j 7 j 7 

This reformulates the van Vleck determinant in terms of a path ordered expo- 
nential 

A 7 (£,y) = s 7 (x,y) d det 



Pexp{- / (V ® t)ds} 



(46) 



This observation serves to illustrate formal similarities (and differences) between 
the van Vleck determinant and the Wilson loop variables of gauge theories. 
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By taking a double derivative with respect to arc length the matrix A(s) 
may be shown to satisfy the second order differential equation 

^A^(s) = - (i^ft") A%. (47) 

See Hawking and Ellis ]rij |, equation (4.21) on page 83, particularized to a 
geodesic flow. The boundary condition on A(s) is that A^ u (s) — » s<5% as 
s — > 0, this being the flat space result for a geodesic spray. This is the promised 
tidal formulation for the van Vleck determinant. 

A particularly pleasant feature of the tidal reformulation is that the case 
of null geodesies can be handled without too much special case fiddling. The 
analogue of the tidal equation is 

^A^(s) = - (R» aa p l a f) AV (48) 

See Hawking and Ellis fl7|| , equation (4.33) on page 88. Note that two of the 
eigenvectors of A suffer trivial evolution. This is a consequence of the two 
normalization conditions on the null tangent vector. Thus A can be thought 
of as either a (d + 1) x (d + 1) matrix with two trivial entries, or as a reduced 
(d— 1) x (d— 1) matrix. The formulations are related by: Au+i) = A(d-i) ©C-fe- 
Other results can be simply transcribed as needed. 

2.7 Reformulation: Formal solution 

To solve the above differential equation (in a formal sense), introduce the one- 
dimensional retarded Green function 

G R (s f ,Si) = {sf-s.jQisf - s^. (49) 

Here O(s) denotes the Heavyside step function. For notational convenience 
define 

Q%{s) = - {R» aiJ pt a tP) . (50) 

Integration of the second-order "tidal" equation for A 7 , with attention to 
the imposed boundary condition, leads to the integral equation 

A%(s) = 8 6" v + [ G B {s lS ')Q\{s')A%{s')ds'. (51) 
Jo 

More formally, one may suppress the explicit integration by regarding Gr(s, s'), 
multiplication by Q(s), and multiplication by s, as functional operators. Then 

(I - G R Q)A = {si}. (52) 
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This has the formal solution 



A = {I-G R Q)- l {aI} 

= (I+[GrQ} + [G r Q] 2 + ---){sI}. (53) 

This formal solution may also be derived directly from equation ( [5l| ) by contin- 
ued iteration. To understand what these symbols mean, note that 

[G R Q]{sI}= f Gnis^'Wvis'Ks'ds' = f\a - s')[Q» v ( s ')]s'ds'. (54) 
Jo Jo 

Similarly 

[G R Q] 2 {sI] = f ds 1 f ds"(s - s')[<rMKs' - s")[Q a v {s")]s". (55) 

Jo Jo 

The formal solution in terms of continued iteration is particularly advanta- 
geous in the case where matter is concentrated on thin shells ]20|, ^lj |2^]. In 
that case the matrix [Q tl v {s)\ is described by a series of delta functions. Due to 
the presence of the Heavyside step function, the formal expansion terminates in 
a finite number of steps 

A = (I + [G R Q] + [G R Q] 2 + ■■■ + [G R Q] N ) {si}. (56) 

Here N denotes the total number of shell-crossings. This type of calculation 
will be explored in considerable detail in a subsequent publication. For the time 
being, one may observe from the above, that when the Riemann curvature is 
concentrated on thin shells, the matrix A is a piecewise linear continuous matrix 
function of arc length. Consequently the reciprocal of the van Vleck determinant 
A 7 (s) _1 = det{A/s} is piecewise a Laurent polynomial in arc length. 

2.8 Weak field approximation 

These formal manipulations permit the derivation of a very nice weak field 
approximation to the van Vleck determinant. Observe 

A 7 (s)- X = det(s- 1 A), 

= expixhi(s- 1 (I-G R Q)- 1 s), 

= exptrln (/ + s'^GrQ^ + s'^GrQ^s + •• ■) , 

= exptr(s- 1 [G fl Q]s + 0(Q 2 )) (57) 

Now, recall that the determinant and trace in these formulae are to be taken in 
the sense of (d+1) x (d + 1) matrices, so in terms of the Ricci tensor 

tr[Q] = - (R aP t a t0) (58) 
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The tr[0(Q 2 )] terms involve messy contractions depending quadratically on the 
Riemann tensor. In a weak field approximation one may neglect these higher 
order terms and write 

A 1 (s)=exp(- S - 1 (G R tr[Q])s + 0{Q 2 )) (59) 

More explicitly 

A 7 (x, y) = exp f ~ J (s - a') (R^rt 13 ) s'ds' + 0( [Riemann] 2 )^ (60) 

This approximation, though valid only for weak fields, has a very nice physical 
interpretation. The implications of constraints such as the (timelike/ null/ 
spacelike) convergence conditions can now be read off by inspection. 

Another nice feature of the weak-field approximation is that the results for 
null geodesies can also be read off by inspection 



A 7 (z,y) = exp ^ J 



C (C - CO {R a pl a l P ) C'dC + 0( [Riemann] 2 ) ) (61) 



Consider, for example, the effect of a thin perfect fluid in an almost flat back- 
ground geometry. The Ricci tensor is related to the density and pressure via 
the Einstein equations, with the result that for a null geodesic 



A 7 (x, y) ps exp ^8ttG J [p+p] 



(C p C ' d C) (62) 



Unsurprisingly, one sees that energy density and/or pressure located at the 
half-way point of the null geodesic is most effective in terms of focussing. 



2.9 Short distance approximation 

Another result readily derivable from the formal solution ( p3[ ) is a short distance 
approximation. Assume that x and y are close to each other, so that s 7 (x,y) 
is small. Assume that the Riemann tensor does not fluctuate wildly along the 
geodesic 7. Then one may approximate the Ricci tensor by a constant, and 
explicitly perform the integration over arc length, keeping only the lowest order 
term in s 7 (x, y) 

A 7 (x,y) =1+1 (R a pt a tP) Sl (x,y) 2 + 0(s 7 (x,y) 3 ). (63) 

This result can also be derived via a veritable orgy of index gymnastics and 
point-splitting techniques Q], see equation (1.76) on page 233. 
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By adopting Gaussian normal coordinates at x one may write 
s~t{x, y) t^{x; 7; x <- y) = (x - yY + 0[s 2 } to yield 

A 7 (x, v) = l + \ {R a /3(x - y) a (x - yf) + 0( Sl (x, y) 3 ). (64) 

In this form the result is blatantly applicable to null geodesies without further 
ado. 

2.10 Asymptotics 

To proceed beyond the weak field approximation, consider an arbitrarily strong 
gravitational source in an asymptotically flat spacetime. Let the point y lie 
anywhere in the spacetime, possibly deep within the strongly gravitating region. 
Consider an otherwise arbitrary spacelike or timelike geodesic that reaches and 
remains in the asymptotically flat region. 

Far from the source the metric is approximately flat and the Riemann tensor 
is of order 0[M/r 3 ]. Using the tidal evolution equation, a double integration 
with respect to arc length provides the estimate 

A"„(«) = (4>)"„ + s{B Q y v + 0[M/s\. (65) 

This estimate being valid for large values of s, where one is in the asymptotically 
flat region. Here (Ao) and (Bo) are constants that effectively summarize gross 
features of the otherwise messy strongly interacting region. For the van Vleck 
determinant A 7 (s) _1 = dct{A/s} = det{B + {A /s) + 0[l/s 2 }}. Thus the van 
Vleck determinant approaches a finite limit at spatial or temporal infinity, with 
A 7 (oo) = det{-Bo} -1 . Finally, define J = (B )- 1 (A ), to obtain 

A ^ (S) = dct{/ + (J/Vo[l/* 2 ]} - (66) 
Which is our desired asymptotic estimate of the van Vleck determinant. 



3 DISCUSSION 

This paper has presented a number of formal developments with regard to eval- 
uating the van Vleck determinant. The evolution of the van Vleck determinant 
as one moves along a geodesic has been studied in some detail. By utilizing 
the Raychaudhuri equation stringent limits have been placed on this evolution. 
There is some hope that it might be possible to formulate more general the- 
orems requiring weaker convergence hypotheses than those discussed in this 
paper. The evolution of the van Vleck determinant has also been reformulated 
in terms of tidal effects. Such a presentation has several technical advantages. 
Amoung other things, it makes manifest the influence of the full Riemann tensor. 
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A formal solution to the resulting integral equation was presented. This formal 
solution in terms of an iterated integral equation has the very nice property of 
terminating in a finite number of steps if the curvature is confined to thin shells. 
Using these formal techniques, a weak field approximation was also developed. 
This weak field approximation depends only on the Ricci tensor. In a similar 
vein, a short distance approximation was developed. Finally, the asymptotic 
behaviour at infinity was investigated. The mathematical machinery developed 
in this paper has a wide arena of applicability. Accordingly, some effort has 
been made to keep the discussion as general as reasonably possible. 
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